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We have calculated the electronic structure of Eu for the bcc, hep, and fee crystal structures for volumes near 
equilibrium up to a calculated 90 GPa pressure using the augmented-plane wave method in the local-density 
approximation. The frozen-core approximation was used with a semi-empirical shift of the f-states energies in 
the radial Schrodinger equation to move the occupied 4f valence states below the F] energy and into the core. 
This shift of the highly localized f-states yields the correct europium phase ordering with lattice parameters 
and bulk moduli in good agreement with experimental data. The calculated superconductivity properties under 
pressure for the bcc and hep structures are also found to agree with and follow a T^ trend similar to recent 
measurement by Debessai et ali 
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A recent paper by Debessai et al. — reports on membrane- 
driven diamond-anvil experiments on europium showing the 
onset and linear increase of superconductivity under hydro- 
static pressure. A transition temperature near 1.8 K was ob- 
served near 80 GPa, increasing to 2.75 at 142 GPa. Previous 
x-ray diffraction experiments by Takemura and Syassen- note 
the bcc, hep, hep-like crystal structures, with transitions oc- 
curring at 12.5 and 18 GPa, respectively. 

Several theoretical works dating back to the mid 1960s 
have focused on the electronic structure of the equilibrium 
bcc structure. Freeman and Dimmock reported the first 
band structure calculations using the non-self-consistent aug- 
mented plane-waves (APW) method.^ Andersen and Loucks 
soon followed with relativistic APW (RAPW) method which 
also was not self consistent and used the same exchange po- 
tential parameter a - \ .0."* Kobayasi et al. performed similar, 
non-self-consistent calculations with the KKR method with a 
- 1.0, 0.9, 0.8 and 0.67.^ Matsumoto et al. related the shape 
of the Fermi surfaces to helical spin ordering then reported the 
first self-consistent KKR results for a = 1.0 and 0.67, placing 
the occupied 4f bands near those from experimenti^. Mat- 
sumoto et al. then reported self-consistent paramagnetic band 
structure using the Xa-KKR method.^ Later Min et al. per- 
formed self-consistent semi-relativistic and fully-relativistic 
linearized muffin-tin orbital (LMTO) total energy calculations 
within the local density approximation (LDA) for ambient 
and high pressure bcc Europium.^ Recently Turek et al. cal- 
culated the exchange parameters, corresponding to magnon 
spectrum, and Neel temperature from a real-space perturba- 
tion approach derived from the tight-binding linear muffin-tin 
(TB-LMTO) method.'" And Kunes and Laskowski calculated 
the magnetic ground state and Fermi surface of bcc Europium 
using the linearized augmented-plane-waves (LAPW) method 
with LDA+U to determine spin spiral states." 

This work reports the results of paramagnetic self- 
consistent, scalar-relativistic, frozen-core augmented plane- 
wave (APW) calculations of the total energies and electronic 
structures of 6s^ valence bcc, fee, and hep europium for 
the Hedin-Lundqvist'^ local-density approximation (LDA) 
method with a semi-empirical shift of the f-state energies. 



In the spirit of previous bcc europium spin-polarized calcu- 
lations that split localized f-states near e/r to determine struc- 
tural and magnetic properties in bcc europium (see for exam- 
ple Refs. i9l10.) , self-consistent scalar-relativistic frozen core 
APW total energy calculations were conducted for the bcc, 
fee, and hep structures with a volume-based empirical shift 
of the f-states to determine the ground state and superconduc- 
tivity properties as functions of pressure. The Brillouin-zone 
summations over the valence states used 285 bcc, 45 hep, and 
505 fee k-points in the iiTeducible portion of the Brillouin- 
zone and a broadening temperature equivalent to 2 mRy. To 
simplify the volume-based shift procedure, the hep structure 
was only considered as ideal with c/a ratio -\/|. Experiment 
however shows a decrease in the c/a ratio from 1.60 at 12.5 
GPa to near 1.56 at 18 GPa.^ 
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FIG. 1 : Europium bcc shifted (solid-line, ■) and unshifted 
(dashed-line, •) f-states total energy showing no minimum. 

In a previous paper on lanthanum Ref. [ij we were able 
to accurately determine ground state and superconductivity 
properties because the unoccupied 4f -bands are localized well 
above the Fermi level. Calculating similar properties for eu- 
ropium is problematic because common approximations to 
density functional theory (DFT) such as the local density ap- 
proximation (LDA) as shown in Fig.[T]fail to reliably capture 
ground state properties due to the highly localized 4f-states 
below the Fermi level. "^ Figure[T]also shows the bcc LDA to- 



tal energies in which we have apphed an empirical shift to 
the f-states as described in detail below. No energy mini- 
mum was found for the unshifted f-states calculation which 
decrease monotonically in energy with volume. 

Since we are most concerned with paramagnetic bulk prop- 
erties and a reliable ground state, the occupied 4f-states near 
the Fermi level were placed below the Fj state and treated as 
part of the core. The remaining two 6p valence electrons per 
atom were treated as scalar-relativistic bands. The volume- 
based shift energy of the f-states was determined in the bcc, 
fee, and hep structures at volume voli with. 



C, 



Ystal 



-'f— shift, vol=i 



((^y/22/3'/') 



(1) 



^ vol, 



where Cc-ysiai was the shift coefficient for the bcc, fee or hep 
structure, voZ, was the volume for the shifted calculation, and 
volref was the reference volume. All structures used the same 
reference volume of 319.20 au^ which is the volume for the 
experimental bcc lattice a = 8.61 au. At this volume the 
bcc shift coefficient Chcc was 0.6897 Ry while the fee and 
hep shift coefficients, C/„ and C^cp respectively, were both 
{2l3fl'^Cbcc which was equivalent to 0.5632 Ry. The 4f- 
energies were shifted in the frozen-core self-consistent iter- 
ations by adding the shift energy {Ef^sHfi, voi=i) to the f-state 
energies on the right side of the radial Schrodinger equation: 



1 d odui 1(1 + 1) 
---(r^-^) + {^^ + V„)u, = Eu,, (2) 

where for I -!>,£- E + E /shift, voi=i- The bcc shifts yield 
the energy minimum as shown in Fig.[T|as the solid-line with 
black squares. 
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FIG. 2: Europium hep (solid-line, ■), bcc (dashed-line, •) 

and fee (dash-dot-line, a) total energies relative to 

fourth-order Birch fit to bcc structure. 

In Fig. 12] we show the total energies for europium hep, bcc, 
and fee using the prescribed shifted f-states as a function of 
volume. The symbols represent the calculated values and the 
lines are from the 4''' order Birch fit. The bcc ground state 



is correctly determined with an equilibrium lattice constant of 
8.25 bohr which is about 4 % below experiment (typical of 
LDA calculations) and a calculated bulk modulus of 14.8 GPa 
that is about 1 % greater than experiment. '"'•'^ The bcc and 
fee total energies are very close with bcc equilibrium about 
0.3 mRy below the fee. 
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FIG. 3: Europium hep (solid-line, ■) and bcc (dashed-line, 
•) enthalpies. 

Since our DFT calculations use temperature T =0 K, the 
enthalpy is equal to the Gibbs free energy G=E -i-pV -TS such 
that 



H(P) = E{P) + PV , 



(3) 



The calculated enthalpies are shown in Fig. [3] as functions of 
pressure for the bcc (circles) and hep (squares) structures. The 
transition pressure occurs where the H(P) curves for the bcc 
and hep cross near 10.6 GPa which is 15.2 % below the transi- 
tion observerd at 12.5 GPa in experiment.^ Figure|4]compares 
the calculated volume ratios of the hep, bcc, and fee as func- 
tions of pressure with the ratios from experiment.'^ Though 
the calculated volume ratios as a function of pressure follow a 
trend similar to experiment they are slightly larger. 

Our bcc results compare well with previous linearized 
muffin-tin orbital, calculations of Min et al!^ Their calcula- 
tions with the paramagnetic 4f electrons as valence states sev- 
erly overestimate the bonding of the localized 4f band, pro- 
ducing a small equilibrium lattice constant 7.61 au and a large 
bulk modulus 45 GPa. Their ferromagnetic calculations split 
the 4f states and generated a lattice constant of 8.24 au and 
bulk modulus of 25 GPa. Placing the 4f-electrons in the core 
significantly improved the bcc equilibrium lattice a - 8.73 au, 
which is about 1.4% above the experimental lattice of 8 .6 1 au 
for both their paramagnetic and ferromagnetic calculations. 
However their paramagnetic and ferromagnetic bulk moduli, 
which are 20 GPa and 19 GPa respectively, are about 33 % 
larger than the 15 GPa found in experiment. Our empirical- 
shift model achieves bcc lattice constants similar to those from 
the ferromagnetic calculations by Min et. al but we also find 
a bcc equilibrium bulk modulus of 14.8 GPa that is very close 
to experiment. 
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FIG. 4: Europium hep (solid-line, ■), bee (dashed-line, •) 

and fee (dash-dot-line, A) volume ratios as funetions of 

pressure. Experimental data (solid-line, D) are from 

Grosshans et al. (Reflla). 
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FIG. 5: Europium shifted f-states bee equilibruim (a=8.25 
a.u.) total and deeomposed density of states. 



The results from our first-prineiples ealculations beeame in- 
put to the rigid-muffin-tin approximation (RMTA) theory de- 
veloped by Gaspari and Gyorffy to yield the high-pressure su- 
pereonducting behavior of bee and hep europium, providing 
additional insight into the high pressure T^ trend of Eu-IIIji2ii^ 



III. DENSITY OF STATES 

The density of states, N(e), were ealeulated by the tetra- 
hedron method'^ using the eonverged self-eonsistent frozen- 
eore bands. The bee and fee N(e) were determined using 285 
and 505 fee k-points respectively. The hep N(e) was cal- 
culated from Gillan's method— as described in Ref. [13| but 
with a larger Fermi temperature of 0.02 Ry to provide addi- 
tional weight to the broadening function within the k-point 
mesh. The sum of the product of energy eigenvalues and ma- 
trix elements using 45 k-points in the irreducible Brillouin 
zone yielded the weighted hep density of states. Figures |5] 
and|6]show the calculated bee density of states for the shifted 
and unshifted 4f-energies respectively at equilibium (a - 8.25 
a.u.). 

The DOS that corresponds to the case of the shifted f-states 
has the Fermi level, ep, located to the right of a peak near 
0.3 Ry and has a large d-like contribution. The DOS of the 
unshifted f-states for a narrow region near ef are shown in 
Fig. |6] The occupied 4f-states are located below e/r and the 
unoccupied 4f-states found above ep. The unshifted f-states 
total DOS at ep, which is composed almost entirely of f-like 
states, is about 100 times greater than the shifted N(ep). 

The significant influence of the d-like states at ef is viewed 
best in Fig. |2l where the ratio N(ef)i/N(ef) is displayed as 
a function of pressure for the hep (a), bee (b), and fee (c) 
structures. The dominant d-like states remain fairly constant 
showing only moderate changes with pressure. The hep d-like 
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FIG. 6: Europium unshifted f-states bee equilibrium (a=8.25 

a.u.) total and decomposed density of states at a narrow range 

near the Fermi level. 



states at ep increase to almost 0.8 near 24 GPa, drop to a low 
of about 0.61 near 70 GPa then increase at higher pressures. 
The bee d-like states at ep decrease from 0.73 at equilibrium 
to 0.68 at near 12.5 GPa. Between 12.5 GPa and 58 GPa the 
bee d-like states increase back to 0.73 then return to about 
0.68 near 58 GPa. At higher pressures the bee d-like states 
increase reaching a maximum of almost 0.78 near 90 GPa. 
And the fee d-like states at ep decrease almost monatonically 
from equilibrium pressure to about 0.60 near 70 GPa. The fee 
d-like states show only minor increases at higher pressures. 

The hep f-like states at ep increase monatonically with pres- 
sure with relatively no change between 25 and 35 GPa and a 
slight decrease after about 85 GPa. At equilibrium pressure 
the hep f-like states are about 0.05 and reach a maximum of 
almost 0.2 near 85 GPa. The bee f-like states at ep increases 
from about 0.05 at equilibrium to almost 0. 1 near 58 GPa, de- 



0.8 



0.6 



0.4 



0.2 



0.8 



0.6 



0.4 



0.2 



s — ^ 

p — X- 

d --^- 

f B- 



~ * - - -X- - - - X ----«- - 




10 20 30 40 50 60 70 80 90 

Pressure(GPa) 

(a) hep ^/(ef )/iV,(ef ) 



' 










s 




















p 


X — 














d 


---.X--- 




0.8 


- 








f 


B 


'X 






-^- 


-X-- 


-^''' 


_^. .---*"-- 






0.6 


- 












- 


0.4 


- 












- 


0.2 


- 










1 


-§ 




^■■B --B B B-B-a-e--Q- 

1 1 1 1 


--X- 

■B- 


:-a-^- 


::»,, 


1 



10 20 30 40 50 60 70 80 90 

Pressure(GPa) 

(b) bcc Ni{£f)/N,(£f) 
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FIG. 7: Europium (a) hep, (b) bcc, and (c) fee angular 
momentum deeomposed DOS divided by total DOS at ef. 



monatonieally with pressure from about 0.035 near 4 GPa to 
almost 0.17 near 90 GPa. 

The hep s-like density of states at £f oscillate with pres- 
sure from a high of almost 0.16 to a low near 0.1. The bcc 
s-like states at £f show only minor changes from equilibrium 
pressure to about 62 GPa. Afterwards they decrease to almost 
0.06 near 90 GPa. The fee s-like states increase monatonieally 
from almost 0.06 near 4 GPa to about 0. 15 slightly with pres- 
sure to equilibrium where is remains relatively constant with 
pressure. 

The hep p-like N(ef) ratio decreases from about 0.18 at 
equilibrium to about 0.09 at 20 GPa where it shows no change 
with pressure. The bcc p-like N{ef) decreases slowly from 
near 0.12 at equilibrium to about 0.06 near 62 GPa, remain- 
ing constant with pressure. And the fee p-like N(e) remains 
constant with pressure from about 0.08 near 4 GPa. 

Overall the Ni{ef)/N,{€f) ratios vary slowly with pressure. 



A. Band Structure 

Figure [8] shows the band structures of bcc equilibrium re- 
spectively with the shift of the energy of the 4f-states (Fig. [Sab 
and without the shift (Fig.[8b]i. The shifted 4f-states are below 
the Fi energy and highly localized within a energy block less 
than 0. 13 Ry. The unshifted 4f-states are still highly localized 
but entangle themselves near the Fermi level, making it im- 
possible to determine the minimum energy for the structure. 

Our empirical-shift model explicitly omits hybridization of 
the occupied 4f levels by placing them below the Fi level, they 
contribute little to chemical bonding.^' Spin-polarized calcu- 
lations by Kunes and Laskowski, Min etai, and others that 
segregate the occupied and unoccupied 4f states respectively 
below and above ef report negligible mixing of the occupied 
4f states.^'" Recent photoemission studies by Wang et al. re- 
ports the highly localized 4f core level at about 0.154 Ry be- 
low the Fermi level. Replicating the 4f core location with such 
proximity to ef would be difficult even with the more sophis- 
ticated LDA+U or the self-interaction corrected (SIC) LSDA 
approaches. 



creases to about 0.07 above 60 GPa, and increases back to 
almost 0. 1 near 90 GPa. The fee f-like states at e/r increase 
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FIG. 8: Europium equilibrium bcc shifted and unshifted f-states band structures. 



IV. SUPERCONDUCTIVITY 

From our densities of states results at the Fermi level we de- 
termine, from first-principles, the electron-phonon coupling 
constant A and the superconducting transition temperature 
Tc- Using McMillan's strong coupling theory,'^ we take the 
electron-phonon coupling constant by the following formula: 



n 



where M is the atomic mass, r\ is the Hopfield parameter— cal- 
culated with the RMTAr^ N(ef ) is the total density of states 
at the Fermi level, ep, < /^ > is the square of the electron-ion 
matrix element at e^ , and < (J- >is the square of the average 
phonon frequency. The Hopfield parameter is determined by 
the quantity 



where the shift energy {Ef^sMft, voi=i) was added to the Fermi 
energy of the 4''' component 1=3 in both Eqns.|2]and[8] Equa- 
tion |6] is used with touching MX spheres to minimize errors 
due to excess interstitial volume as demonstrated in Ref. |2g, 
I27I and|T3l The APW total energy results were used to calcu- 
late the pressure variation of the bulk moduli B. These values 
were used to determine the pressure variation of the Debye 
temperature 6d, which has no experimentally measured bcc 
equilibrium 6d below a temperature of 100 K,^-' using the for- 
mula by Moruzzi et al.'J^ 



213.4V— 



(9) 



where B is the bulk modulus (in GPa), rg is the Wigner-Seitz 
radius in Bohr units, and M is the atomic mass. The approxi- 
mate average phonon frequency < w^ > is determined as 



77 = N(ef) <I>, (5) 

using the RMTA. The square matrix element < fi >, which is 
derived from multiple-scattering theory and is also determined 
by the RMTA proposed by Gaspari and Gyorffy, is defined by 
the formula 



< a? >= ]^el , (10) 

and is used both in the denominator of Eq.|4]and in the pref- 
actor of the McMillan equation for the determination of the 
superconductivity transision temperature T^^ 



<f >= 



fF 



2(l+l)sin\6M-6,)N,NM 






(6) 



where 5i are scattering phase shifts, Ni(ef) is the /''' compo- 
nent of the DOS per spin, and A^, ^ is the free-scatter DOS, 

/ r-R, 

N\^^ ^ J-L(2l + \) \ uhr,efydr. (7) 

7!" Jo 

The scattering phase shifts 61 were calculated from the ra- 
dial wave functions «/, spherical Bessel functions ji and Neu- 
mann functions «/ at the muffin-tin radius /?, from the formula 



u'i(r,eF) _ j'ljkr) - n'i(kr)tan 61 
ui(r, £f) jiikr) - ni(kr)tan 5i 



(8) 



< ClJ> 



-exp[- 



-1.04(1 + /i) 



(11) 



1.2 '^^-yu*(l-H0.62i) 
where the Coulomb pseudopotential yU* was fixed at 0.13 for 
all three structures. 

From our total energy and density of states calculations we 
determined the input parameters for the Tc calculations for 
the shifted 4f-states for the bcc, hep, and fee structures. The 
bulk modulus for Eu bcc is shown in Figure |9a| as a function 
of pressure. The average phonon frequencies show a similar 
trend, as displayed in Figure |9bl since they are derived from 
the bulk moduli. These trends, where B is directly related 
to < w >, also agree qualitatively with the proportionality 
relations used for the lighter elements of fee yttrium,^'' fee 
scandium,^^ and lanthanum^ 
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FIG. 9: Bulk modulus B and average phonon frequency < w > of hep (solid-line, ■), bcc (dashed-line, •) and fee 

(dash-dot-line, k) as functions of pressure. 
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FIG. 10: Total DOS per spin at e^ and electron-ion matrix element < fi > of hep (solid-line, ■), bcc (dashed-line, •) and fee 

(dash-dot-line, a) as functions of pressure. 
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6i) for / - (squares), / = 1 (circles), and / = 2 (triangles) and the total and 
decomposed Hopfield parameter 77 of bcc europium as functions of pressure. 



The density of states at the Fermi level N(ep), electron-ion 
matrix element < /^ >, bcc phase shift factor, and bcc total 



and decomposed Hopfield parameter rj as functions of pres- 
sure are shown in Figs. llOal 11 Obi lllal and lllbl respectively. 
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FIG. 12: Hopfield parameter r\ and electron-phonon coupling constant X of hep (solid-line, ■), bcc (dashed-line, •) and fee 

(dash-dot-line, a) as functions of pressure. 



After about 20 GPa the density of states varies slowly with 
pressure to 90 GPa for the hep and fee. The bcc N(ef^) in- 
crease rapidly after about 55 GPa from near 3.7 states/Ry/spin 
to almost 5.8 states/Ry/spin near 70 GPa. The hep and fee 
N(ef ) also decrease after 70 GPa. 

Although the N(ef^) vary slowly with pressure, the electron- 
ion matrix element < /^ > in Eqn. |5] increases rapidly with 
pressure in all three structures. The increase of the electron- 
ion matrix element < fi > comes primarily from the large 
p-d phase shift factor which is near one at high pressures as 
shown in Fig. I Hal for bcc. Figure lllbl shows the bcc total 
and decomposed Hopfield parameter as a function of pressure. 
The large increase in total // with pressure comes from the 
large p-d 77 contribution. The hep and fee also have large p-d 
phase factor and 77 components that increase the total 77. 

Figures |12a| and |12b| respectively show the total 77 and 
electron-phonon coupling constant A as function of pressure 
for the hep, bcc, and fee. The increasing 77 values in the nu- 
merator of Eqn. ID along with moderate phonon frequencies 
in the denominator help explain the increase of the electron- 
phonon coupling constant A with pressure. Above 55 GPa the 
bcc and hep A goes above 0.4 producing an increase in the 
superconductity transistion temperature Tc shown in Fig. [T3] 
The calculated bcc and hep Tc flattens near 80 GPa because 
the A levels off at about 0.53 near the same pressure range. 
The fee shows a peak in Tc near 70 GPa but decreases at higher 
pressures due to the decreasing fee 77 that yields smaller A with 
pressure. 
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V. CONCLUSIONS 

In conclusion, the electronic structure and superconduc- 
tivity properties of europium have been studied as functions 
of pressure using the Hedin-Lundqvist LDA form of the 
exchange-correlation and shifted 4f-band energies within the 
APW method. In particular, we investigated the total energies 
for the bcc, fee, and hep structures; the density of states, for 
the bcc and hep structures and the band structure for the bcc 
and hep; and the superconducting properties for the bcc and 
hep structures through high pressures, using the RMTA. 

On the basis of our total energy calculations, the bcc lat- 
tice has been shown to be the equilibrium structure, while the 



hep becomes favorable at high pressures. Our density of states 
calculations agree well with previous results for the bcc struc- 
ture. Our calculations of the superconducting properties show 
a gradual increase of the electron-phonon coupling constant to 
values around 0.5, reproducing the onset of superconductivity 
at pressures between 60 and 90 GPa. 
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